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ABSTRACT 


This  report  is  a  theoretical  study  of  the  errors  due  to  ran¬ 
dom  noise  to  be  expected  in  the  measurements  of  satellite  posi¬ 
tion,  velocity,  and  other  orbital  parameters  obtained  by  the  U.S. 
Naval  Space  Surveillance  System.  Both  phase  and  phase-rate  data 
are  considered  with  a  view  toward  the  following  applications: 

(a)  weighting  data  from  satellite  passes  received  at  several  sta¬ 
tions  to  produce  better  estimates  of  both  position  and  velocity, 

(b)  improving  first  pass  capability  to  determine  orbital  parame¬ 
ters,  and  (c)  providing  real-time  designation  data  for  tracking 
systems.  Results  are  displayed  as  equal  error  contours  on  ver¬ 
tical  plots  of  the  Spasur  fence  for  basic  position  and  velocity 
components. 


PROBLEM  STATUS 

This  is  a  final  report  on  one  phase  of  the  problem;  work  on 
other  phases  continues. 


AUTHORIZATION 

NRL  Problem  R02-35 
Project  RT AD-22-001/652- 1/5434-00-000 


Manuscript  submitted  October  E6,  1962. 
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A  THEORETICAL  ERROR  ANALYSIS  OF  THE  U.S.  NAVAL 
SPACE  SURVEILLANCE  SYSTEM 


INTRODUCTION  AND  DEFINITION  OF  TERMS 
Description  of  the  Spasur  System 

The  U.S.  Naval  Space  Surveillance  (Spasur)  System  is  a  multi-static  cw  radar  sys¬ 
tem  designed  for  the  detection  and  identification  of  satellites.  A  satellite  is  identified 
principally  by  the  values  of  its  orbital  parameters.  At  present  an  observation  is  said  to 
be  identified  when  it  correlates  within  predetermined  tolerances  to  predictions  based  on 
previously  determined  orbital  elements.  For  a  new  satellite,  orbital  elements  are  now 
determined  from  the  positions  on  consecutive  passes.  However,  the  system  is  capable  of 
providing  limited  velocity  information  which  may  be  used  on  the  first  pass  of  a  satellite 
to  compute  orbital  elements.  Measurements  of  position  and  velocity  contain  errors, 
which  are  themselves  functions  of  satellite  position  and  velocity,  so  that  the  Spasur  Sys¬ 
tem  is  more  accurate  in  some  regions  in  its  beam  than  in  others.  In  this  first  main  sec¬ 
tion  we  shall  describe  how  these  errors  arise,  in  the  second  section  we  derive  formulas 
which  express  these  errors  as  functions  of  satellite  position  and  velocity,  and  in  the  third 
section  we  present,  as  equal  error  contours,  the  results  of  some  computations  showing 
the  variation  in  error  with  the  position  in  the  beam. 

For  our  purposes  the  following  simplified  model  of  the  system  will  be  sufficient.  The 
system  consists  of  three  cw  transmitters  and  four  receivers  arranged  along  a  great  cir¬ 
cle  inclined  at  33°57'  to  the  equator  (Fig.  1).  The  receivers  are  labeled  1,  2,  3,  and  4 
from  west  to  east.  The  plane  of  the  great  circle  is  called  the  "fence  plane,"  or  simply 
"fence."  Each  receiver  consists  of  one  interferometerlike  phase  measuring  device 
oriented  in  the  fence  plane  and  another 
oriented  perpendicular  to  it  (1).  When 
a  satellite  passes  through  the  fence, 
its  zenith  angle,  or  "space  angle,"  d 
at  each  receiver  station  is  deter¬ 
mined  by 

sin  ^  T  (1) 

where  ^  is  the  radar  wavelength,  d  is 
the  baseline  or  spacing  between  the 
antennas  of  the  interferometer,  and 
flp  is  the  phase  difference  between  the 
signals  received  by  the  antennas  of 
the  interferometer  (in  rotations).  It 
can  be  seen  that  a  two- station  fix  will 
determine  the  position  of  a  satellite 
(by  triangulation).  Also,  because  of 
errors  in  the  measured  values  of 
a  three- station  fix  will  define  three 
positions  (Pj2(  P23»  and  Pjj  in  Fig. 

1),  while  a  four-station  fix  will  define 
six  positions. 


/ 

/  / 
// 
// 


/A”” 
/  /  ' 


Fig.  1  -  Spasur  fence  plane,  showing 
a  three-station  fix 
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The  transmitters  have  a  narrow  beamwidth  perpendicular  to  the  fence,  so  that  a  sat¬ 
ellite  is  observed  for  a  finite  period  of  time.  This  enables  us  to  estimate  phase  rates 
and  therefore  angle  rates  The  angles  and  angle  rates  from  two  stations  can  be  used 
to  estimate  the  component  of  a  satellite's  velocity  in  the  fence  plane.  The  interferome¬ 
ter  oriented  perpendicular  to  the  fence  can  be  used  to  obtain  estimates  of  velocity  com¬ 
ponents  perpendicular  to  the  fence,  and  also  times  of  fence  crossings. 


Notation  for  the  Statistical  Concepts  Used 

In  this  report  either  <ii>  or  n  will  denote  the  expected  value  of  a  random  variable 
u.  Then  <u^>  or  is  the  mean  square  value  of  u,  and  /<u7>  or  is  the  root 

mean  square  (rms)  value  of  u.  When  u  has  zero  mean,  <u2>  and  \/<ii2>  are  frequently 
called  the  variance  and  standard  deviation  of  n .  The  covariance  between  two  random 
variables  u  and  v  is  <uv>,  and  <uv>  =  o  if  u  and  v  are  statistically  independent.* 

Now  suppose  f(ii)  is  a  function  of  »i,  where  u  is  a  quantity  which  is  to  be  experi¬ 
mentally  determined.  Then  a  small  error  in  the  measured  value  of  u  will  produce 
an  error  in  f(ii>  equal  to 


sf  =  f'(u)  Sii  (2) 

where  f'(u)  is  to  be  evaluated  at  the  true  value  of  u.  Equation  (2)  is  valid  if  f '(u)  is 
continuous  at  this  point.  In  this  report  errors  in  the  independent  variables  will  always 
be  random  variables  with  zero  mean,  so  that  the  mean  square  error  and  rms  error  in  f 
are 

<8f2>  =  [f'(u)]*  <i>u*>  (3) 

v/<8f *>  =  |f'(u)l  \/<5u2>  .  (4) 

Similarly,  if  f  =  f(u,v) ,  then  small  errors  5u  and  8v  will  produce  the  error 

*' =  1^ «■> 

and  we  have 

=  (f )'  .  2  |t  f  <*„  Sv>  ,  (|1)'  <iv2>  (}■) 

=  V  (ft)' 

It  should  be  emphasized  that  Su,  6v,  and  8f  are  random  variables,  while  <Su*>,  <Sv2>, 
and  <bf^>  being  variances  are  properties  of  the  distribution  function  and  are  not  random 
variables.  (For  example,  <»f*>  is  in  general  a  function  of  u  and  v  which  does  not  con¬ 
tain  random  elements.) 


Determination  of  Angle  and  Phase  Errors 

A  satellite's  position,  velocity,  and  orbital  parameters  will  be  computed  from  the 
angles  •'  and  angle  rates  0;  therefore  from  Eq.  (1)  any  errors  in  the  information  provided 


'■See  Reis.  2-4  for  an  exposition  of  these  statistical  concepts. 
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<5^^>  =  sec  6^  <hrf^>  .  (®) 

Equation  (5)  is  a  good  first  order  approximation  so  long  as  sec  9  is  not  too  large  (i.e., 
so  long  as  the  satellite  is  not  too  near  the  horizon  of  a  receiver).  This  same  condition 
must,  of  course,  hold  true  for  Eq.  (6)  to  be  valid,  with  the  further  condition  that  &iii  never 
be  much  smaller  than  i<P.  This  last  condition  holds,  since  it  can  be  shown  that  Sq)  is  of 
the  same  order  of  magnitude  as  J<p/T,  where  T  is  the  satellite's  time  of  passage  and  is 
generally  less  than  2  seconds. 

Now  the  phase  <p  is  measured  and  recorded  continuously  by  a  "phasemeter,"  whose 
output  during  the  time  of  a  satellite's  passage  through  the  fence  is  approximately  a 
straight  line  disturbed  by  additive  noise.  We  shall  assume  that  the  phase  <p  and  phase 
rate  <p  will  be  estimated  by  fitting  a  straight  line  to  the  phasemeter  output  by  least 
squares  (or  some  similar  method).  The  quantities  <Sip2>  and  <Sq)2>  can  then  be  esti¬ 
mated  from  the  theory  of  least  squares  if  we  know  the  variance  of  the  noise  disturbing 
function,  but  this  latter  quantity  can  be  estimated  from  the  mean  square  value  of  the 
residuals  (difference  between  actual  values  and  smoothed  values  of  data). 

In  other  words,  once  a  data  smoothing  method  has  been  selected,  the  values  of  the 
parameters  <s<p^>  and  can  be  estimated  from  the  phasemeter  data  alone,  without 

any  additional  Information  about  the  satellite.  This  is  very  important,  since  the  true 
values  of  <Snp^>  and  <si*>  depend  in  a  complicated  way  on  the  signal-to- noise  ratio, 
which  in  turn  depends  upon  the  satellite's  position,  radar  cross  section,  and  the  system 
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noise.”  Consequently,  in  what  follows,  we  shall  consider  and  to  be  inde¬ 

pendent  of  any  other  satellite  parameter. 


Uses  of  the  Errors  to  be  Computed 

The  errors  to  be  computed  are  the  following:  (a)  east-west  position  errors,  (b)  north- 
south  position  errors,  (c)  north-south  velocity  errors,  and  (d)  errors  in  east-west  and 
total  velocity  and  certain  orbital  parameters  for  a  set  of  satellites  in  circular  orbits. 
These  errors  are  defined  in  greater  detail  in  the  next  section.  "East-west"  refers  to 
components  of  position  or  velocity  in  the  fence  plane,  while  "north- south"  refers  to  com¬ 
ponents  perpendicular  to  the  fence.  Some  of  the  uses  of  these  errors  are  the  following: 

1.  Determination  of  errors  in  predicted  future  positions  and  orbit  parameters.  Each 
orbit  parameter  is  a  function  of  six  quantities;  the  three  position  coordinates  and  the 
three  velocity  components  at  the  time  of  fence  crossing.  The  error  in  an  orbit  parame¬ 
ter  can  then  be  computed  from  a  generalization  of  Eq.  (4').  Examples  are  given  later 
(see  page  22). 

2.  East- west  position  estimates  for  a  three- station  or  four- station  fix.  As  shown  in 
Fig.  1,  a  three-station  fix  will  define  three  estimates  of  position,  while  a  four-station  fix 
will  define  six  estimates  of  position.  The  question  arises  as  to  how  to  use  these  estimates 
to  form  a  good  estimate  of  position.  (Also,  the  theory  described  here  is  used  for  north- 
south  position  and  velocity  estimates.) 

For  each  position  coordinate,  say  r,  we  have  the  three  estimates  rjj,  rjj,  rjj 
obtained  from  the  station  pairs  1-2,  2-3,  3-1.  As  an  estimate  for  r  we  might  use  the 
simple  arithmetical  mean  ?  =  (^12*  ^23*  estimate  has  the  disadvantage 

that  it  gives  equal  weight  to  all  the  estimates,  when  in  fact  some  of  these  estimates  may 
be  much  better,  or  worse,  than  others.  If  we  wish  to  give  proper  weights  to  the  different 
estimates,  we  might  use  the  following  general  theorem:  ^ 

If  xj,  xj,  X3  are  three  independent  estimates  of  a  variable  x  with  zero  mean  errors 
and  rms  errors  of  o(x^),  ^(Xj),  jfxj),  then  the  best  linear  estimate  of  x,  i.e.,  the  esti¬ 
mate  of  the  form  \i  xj  +  *  ^3 ’‘a  (where  the  V’s  are  constants)  which  has  zero  mean 

and  minimum  rms  error,  is  given  by 


*1  ^  ^  ^a 

cr2(Xj')  cr^CXj) 
cr2(Xj)  cr^fxj) 


The  rms  error  of  x  is  given  by 


cr(x)  = 


/ 


a^(  Xj) 


CT^fXj) 


(7) 


(8) 


In  order  to  see  how  the  estimate  given  by  Eq.  (7)  compares  with  the  arithmetic  mean, 
we  consider  the  case  for  two  estimates.  We  have 


'^Two  other  parameters  which  affect  the  phase  errors  are  phasemeter  bandwidth  B  and 
duration  of  sij;nal  T.  The  phase  and  phase  rate  errors  vary  inversely  as  BT  and  BT^ 
respectively.  This  dependence  on  T  implies  a  dejjendence  on  satellite  \eiotity  and  height. 
This  is  a  theorem  in  the  theory  of  estimation,  and  can  be  found  in  most  elemimtary  or 
intermediate  statistics  texts. 
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*1 

*2 

T 

B  - 

1 

1 

A  " 

cr^fXj) 

a^(i)  = 

1 

1 

^*(Xi) 

so  that 


<T(x)/arx) 


1 

2  \o^(Xj)  <7(Xj)  j 


(9) 


From  this  last  equation  we  get  <y(x)/a(x)  >  i  (since  q  >  0  implies  q+  i/q  ^  2),  where 
the  equality  sign  holds  only  if  c^fxj)  =  o^fxj). 


Now,  in  our  case,  for  xj,  xj,  xj  in  Eq.  (7)  we  have  rjj,  rjj,  rjj,  and  for  a^(x^), 
ct*(x2),  we  have  the  quantities  computed  by  replacing  f  by  rjj,  T23,  rjj  in  Eq, 

(3'),  The  development  of  ^(r)  in  terms  of  system  input  errors  is  the  subject  of  the  next 
section  (see  especially  Eqs.  (13)). 


Also,  since  the  errors  in  rjj,  rja,  are  not  independent,  the  conditions  for  the 
theorem  cited  above  are  not  fulfilled,  so  that  the  estimate  given  by  Eq.  (7)  is  not  opti¬ 
mum.  The  best  estimate  is  a  more  complicated  expression  containing  crosscorrelation 
terms.  Equation  (9)  still  holds,  however,  and  x  is  always  a  better  estimate  than  x. 


DEHNITION  OF  ERRORS  AND  COMPUTATION  FORMULAS 

In  this  section  we  define  the  errors  described  in  the  preceding  subsection  and  derive 
the  computation  formulas  for  their  calculation.  In  the  next  main  section  we  present  the 
results  of  some  calculations. 


East- West  Errors:  Coordinate  Systems 

At  the  present  time  the  actual  input  position  parameters  to  the  Spasur  data  process¬ 
ing  system  are  earth-centered,  earth-fixed,  polar  coordinates  (r,  i).  For  purposes  of 
calculation  it  was  found  convenient  to  impose  a  rectangular  coordinate  system  for  each 
pair  of  stations  (Fig.  2).  The  x-axis  (see  Fig.  2)  is  the  line  of  symmetry  between  the 
two  stations  (Sj  and  Sj)*  which  are  separated  by  an  earth- centered  angle  of  2a.  The 
y-axis  points  in  a  westerly  direction.  The  sign  convention  is  that  angles  are  measured 
from  the  vertical  and  are  positive  when  measured  counterclockwise.  Thus  in  the  figure 
^1  =  -52  degrees  and  =  '*■72  degrees.  The  radius  of  the  earth  is  a,  and  ^l^  and  ^2 
are  the  slant  ranges  of  the  satellite  from  the  stations. 

Each  point  in  the  fence  plane  is  defined  by  angular  coordinates  (dj,  93),  rectangular 
coordinates  (x,  y),  and  polar  coordinates  (r,  C).  We  have  the  following  transformations: 

=  x*  +  y* 

C  = 


arc  tan  (y/x) 


(10) 


6 


NAVAL  RESEARCH  LABORATORY 
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f  S 

OP  :  J  X*  +  V* 

*  a 

{ ■ 

ore  ten  y/i 

*z- 

Xl- 

a  cot  a 

AX, 

■  x-x, 

Y|‘ 

0  «m  a 

AY, 

.  Y-Y| 

Fig.  2  -  Geometry  of  the  fence  plane 


where 


X 


y 


_ a _ 

tan  -  tan  u'2 

_ a _ 

tan  j  -  tan  c  ^ 


(  sin 

(sin 


sec  -  sin  ^2  '^2'^ 

sec  i/'j  tan  ^2  ~  ^2  '^2 


(11) 


i/'j  -  +  a 

4)2  =  fi'2  -  a  • 

We  can  also  write  these  last  equations  as  follows: 

X  =  — ; - - -  (sin  6^,  cos  i/',  -  sin  cos  41,) 

sin(0i-y,2)  ^  ‘ 


y  =  — : — :: -  (sin  sin  i/*,  -  sin  6/,  sin  4't)  ■ 

sin  (vj  -  V  2)  *  * 


(11') 


Definition  and  Calculation  Formulas  for  East- West  Errors 

The  three  east-west  errors  to  be  calculated  are  the  errors  in  the  input  parameters 
<'r^>  and  <^4^^  and  another  error  <i>w2>  defined  by 
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<dv^>  -  <bx^>  +  <by^>  .  (12) 

Since  the  height  of  a  satellite  is  r  -  a,  <3t^>  may  be  interpreted  as  the  mean  square 
error  in  height.  ■  ')w2>  may  be  interpreted  as  the  expected  value  of  b?  •  Sr,  where  s?  is 
the  position  error  vector  pointing  from  the  true  position  to  the  measured  position.  (We 
have  r  ;  [x.y],  3r  =  [Sx,  e>y] ,  and  s?  •  b?  -  bx^  +  Sy^,)  <bv^>  is  perhaps  the  best  sin¬ 

gle  measure  of  the  positional  accuracy  of  Spasur. 

Applying  ^  to  both  sides  of  Eq.  (10),  we  get  2r  Sr  =  2x  Sx  +  2ySy.  Squaring  and 
applying  <>,  we  get 


<Sr^>  =  — «  (x^  <Sx^>  +  y^  <Sy2>  +  2xy  <Sx  Sy>')  . 


Similarly, 


<bl^>  =  (y^<Sx*>  +  x^  <Sy2>  -  2xy  <Sx  Sy>)  . 

From  Eqs.  (12)  and  (13)  we  get 

<Sw^>  =  <Sx^  +  Sy*>  =  <Sr*>  +  <&i^>  . 


(13a) 


(13b) 


(14) 


Now  all  we  have  to  do  is  to  express  <bx^>,  <Sy*>,  and  <SxSy>  as  functions  of  8^  and 
■>2.  The  derivation  is  given  below,  and  we  get  the  following  results  expressed  in  matrix 
form: 


<ix^> 

1  1 

<ixiy> 

1 

(tan  1^1  -  tan  <p2')^ 

tan  1^2 

sec*  4>i  <b9^> 

sec*  ^>2  <5^2*^ 

1 

A 

CB 

>v 

V 

_ 1 

tan*  1^2  tan* 

(15) 


— 

“  ” 

1 

<Sx*> 

1  1 

1 

cos*  02  <80j^> 

<5x8y> 

sin*  (i/'i  -  i/'2) 

tan  02 

cos*  0j  <hd2> 

<Sy*> 

tan*  02  tan*  0j 

These  equations  can  be  obtained  by  applying  b  to  both  sides  of  Eq.  (11)  multiplying 
the  results,  and  applying  <>.  This  is  rather  laborious,  and  a  simpler  derivation  follows. 
The  equations  of  the  two  lines  of  sight  are  (see  Fig.  2) 


Ayj  -  Axj  tan  =  0,  Ayj  -  Ax2  tan  i/;2  “  0 

with  Axj  =  x-xi,  Ayj^y-yj,  etc.  The  point  coordinates  x  and  y  must  satisfy  these 
equations  identically,  so  that  we  may  apply  &  to  obtain 


by  -  bx  tan  =  Axj  sec^ '/'j  S0j  =  juj  sec  i/<j  bfc'j 


by  -  bx  tan  <^2  -  ^2  “  ^2  '^2  ^^’2' 


Solving  for  bx  and  2>y  we  get 
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Sx  = 


h  = 


tan  -  tan  02 

_ 1 _ 

tan  0j  -  tan  02 


(^2  wc  02  St?2  ■  ^1 


(M2  *86  02  tan  0|  S(?2  ”  Mi  »ec  0i  tan  02  S^i) 


Equation  (15)  can  now  be  obtained  applying  <>  to  the  products  Sx  Sx,  Sx  Sy,  &y  Sy.^  Equa¬ 
tions  (15),  (14),  (13)  and  (5)  are  the  computation  formulas  for  <ir^>,  and  <Sw^>. 


Further  Important  Properties  of  the  East-West  Errors 

In  what  follows  we  assume  that  (X/d)^  <S(p2>  has  the  same  value  at  all  stations,  so 
that,  from  Eqs.  (15),  (14),  (13),  and  (5),  <Sr2>,  <5C*>,  and  <iw*>  are  proportional  to 
(\/d)2  <8(p*>.  We  set  (X/d)^  <8(p^>  =  1 ,  or  equivalently,  =  aec^9,  so  that  the  tabu¬ 
lated  values  of  v/<Sr2>,  and  y/<Sw2>  will  be  equal  to  the  true  values  upon  multipli¬ 

cation  by  the  true  value  of  (X/d) 

Now  suppose  we  make  the  same  assumptions  about  (X/d)^  <Sip^>,  and  set  (X/d)^<S(p^>=  l. 
Then  <Sr*>,  <si^>,  and  <Sw*>  will  be  proportional  to  <8r*>,  <S4*>,  and  <Sx*>  +  <Sy^>, 
assuming  no  error  in  position.  For  example  the  radial  component  of  velocity  is 


SO  that 


^  *^2)  +  pi  S  (^)  +  ^2  S  j- 

If  we  assume  no  error  in  position,  we  have  S(dr/d0^)  =  S(dr/d02)  =  0,  so  that 


But 


Sr 


-  _5l 


S0,  ♦  ^  S0 


(16) 


so  that,  from  Eqs.  (5)  and  (6),  V<8r*>  and  >/<8f*>  differ  only  by  a  multiplicative  constant, 
namely,  y<8<p*>/V<S'i>^>* 

Referring  again  to  Eq.  (16)  we  see  that  <8r^>  is  proportional  to  that  part  of  the  mean 
square  error  in  r  which  is  independent  of  velocity.  The  total  error  in  r  depends  on  four 
independent  terms,  c^,,  Pj,  0,,  and  ffjf  would  be  hard  to  tabulate  as  a  function  of 
position  and  velocity. 

Using  the  same  reasoning  as  in  the  paragraph  above,  and  assuming  no  error  in  posi¬ 
tion  we  see  that  <8w^>  is  proportional  to  the  expected  value  of  8>f  .  8(f,  where  8v  is  the 
velocity  error  vector.  The  quantity  y<8v  •  8v>  is  not  the  expected  error  in  the  length  of 
the  velocity  vector.  We  may  obtain  the  latter  quantity  as  follows.  <8r3>  is  proportional 
to  the  mean  square  error  of  the  velocity  component  in  the  direction  of  r,  and  <SC^>  is 


’■■'Since  and  ^2  measured  at  different  stations,  we  are  assuming  <Sffj  S02>  =  0- 
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proportional  to  the  mean  square  error  in  the  velocity  component  perpendicular  to  r.  Now 
IC?|2  =  +  r2  so 

<S|7|*>  =  -JL  (r2  <Sr2>  +  <8j*>)  . 

171 ' 

This  means  that  the  equal  error  contours  for  position  error  may  be  used  to  determine 
velocity  error,  assuming  no  error  in  position.  Multiplying  a  value  from  Figs.  3,  4,  5  by 
(K/d)  cr(ii),  instead  of  (X^d)  cr(<p),  we  obtain  a  velocity  error  (instead  of  a  position  error). 
Then  using  the  measured  values  of  r ,  f ,  and  4  and  the  above  formula  we  may  determine 
i7(v).  Another  approach  to  this  problem  is  given  later  (see  Eq.  (24)). 


North-South  Errors 

Once  the  location  of  a  satellite  in  the  fence  has  been  determined  from  east-west 
data,  the  north- south  phase  and  phase  rate  from  any  one  station  can  be  used  to  determine 
both  the  velocity  perpendicular  to  the  fence  and  the  north- south  position.  (At  a  given 
time  the  satellite  may  not  be  exactly  in  the  fence,  but  some  distance  away  from  it.)  Errors 
in  these  quantities  are  used  to  compute  errors  in  velocity  and  orbital  elements  and  in  the 
time  of  fence  crossing.  In  this  subsection  we  treat  position  and  time,  and  in  the  next  sub¬ 
section  we  will  consider  the  use  of  phase  rate. 

Let  V  be  the  ^ace  angle  in  the  north-south  direction.  As  before  we  have  Xip=  d  sin  v, 
whence  a(-n)  =  (X/d)  (l/co%  Since  rj  -  o  in  the  fence,  we  have  cos  tj  «  i,  so 

<y(T?)  =  .  (17) 

Now  if  r,  is  the  distance  from  station  i  to  the  satellite  and  ^n^  is  the  angle  meas¬ 
ured  at  that  station,  then  the  corresponding  distance  a^  from  the  satellite  to  the  fence  is 
s,  =  The  error  is  given  by 


and  since  Vi  •  0  we  have 


t7(sj)  =  r^  ff(vi)  ■ 


(18) 


To  combine  data  from  different  stations  we  use  Eq.  (7).  For  two  stations  this  gives 
an  estimate  s  with  variance 


a*(s) 


(19) 


where 


=  (r(r?|)  =  crfTjj) 


Notice  that  we  may  normalize  this  formula  if  we  measure  all  distances  in  units  of 
distance  between  the  two  stations;  i.e.. 


a2(s)  _  (ri/a)^ 

(rj/B)*+ (r^a)* 
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Fig.  3a  -  Equal  error  contours  of  »  10*^:  stations  1  and  2.  The  curve  parameters  may  be  interpreted  as  error 

in  statute  miles  assuming,  for  example,  phase  data  having  18  degrees  (50  counts)  rms  phase  error  from  a  455-ft  (50  K) 
baseline. 
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Fig.  3b  >  Equal  error  contours  of  x  10*^:  stations  2  and  3.  The  curve  parameters  may  be  interpreted  as  error 

in  statute  noiles  asstiining,  for  example,  phase  data  having  18  degrees  (50  counts)  rms  phase  error  from  a  455-ft  (50  K) 
baseline. 


Fig.  3c  -  Equal  error  contoux’e  of  ^<Sr^  x  10*^:  stations  1  and  3.  The  curve  parameters  may  be  interpreted  as  error 
in  statute  miles  assuming,  for  example,  phase  data  having  18  degrees  (50  counts)  rms  phase  error  from  a  455-ft  (50  \) 
baseline. 


Fig.  3d  -  Equal  error  contours  of  ^<8r*>  *  10'*:  stations  1  and  4.  The  curve  parameters  may  be  interpreted  as  error 
in  statute  miles  assuming,  for  example,  phase  data  having  18  degrees  (50  counts)  rms  phase  error  from  a  455-ft  (50  K) 


Fig.  4a  -  Equal  error  contoura  of  ^<Sl^ :  etations  1  and  2.  The  curve  parameters  may  be  interpreted  as  error  in 
earth-centered  angle  measured  in  the  Spasur  plane  in  units  of  thousandths  of  a  degree, assuming  phase  data  having,  for 
example,  18  degrees  (50  counts)  rms  phase  error  from  a  455-ft  (50  \)  baseline. 
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Fig.  4b  -  Equal  error  contoure  of  etatione  2  and  3.  The  curve  parameters  may  be  interpreted  as  error  in 

earth-centered  angle  measured  in  t^  Spasur  plane  in  units  of  thousandths  of  a  degree, assuming  phase  data  having,  for 
example,  18  degrees  (50  counts)  rms  phase  error  from  a  455-ft  (50  K)  baseline. 
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Fig.  4d  -  Equal  error  contours  of  ^<S^*>z  stations  1  and  4.  The  curve  paranneters  may  be  interpreted  as  error  in 
earth-centered  angle  measured  in  the  Spasur  plane  in  units  of  thousandths  of  a  degree, assuming  phase  data  having,  for 
example,  18  degrees  (50  counts)  rms  phase  error  from  a  453-ft  (50  X.)  baseline. 
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5«  -  Equal  error  contoura  of  »  10'*:  etatione  1  and  2.  The  curve  parameters  may  be  interpreted  as  error 

statute  miles  assuming,  for  example,  phase  data  having  18  degrees  (50  counts)  rms  phase  error  from  a  455-ft  (50  K) 
.seline. 
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Fig.  5b  -  Equal  error  contours  of  »  10'^:  stations  2  and  3.  The  curve  parameters  may  be  interpreted  as  error 

in  statute  miles  assiuning,  for  exanq>le,  phase  data  having  18  degrees  (50  counts)  rms  phase  error  from  a  455>ft  (50  K) 
Wseline. 
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Fig.  5d  -  Equal  error  contours  of  ^<8w*>  »  10'^:  stations  1  and  4.  The  curve  parameters  may  be  interpreted  as  error 
in  statute  miles  assiiming,  for  example,  phase  data  having  18  degrees  (50  counts)  rms  phase  error  from  a  455-ft  (50  X.) 
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This  means  that  level  curves  for  the  error  function  at  different  pairs  of  stations  will  look 
the  same,  except  for  a  constant  magnification.  Some  level  curves  are  plotted  in  Fig.  6. 

In  measuring  the  time  t  when  the  satellite  is  at  a  distance  s  there  will  be  an  error. 
We  have  t  =  s/v  from  which  we  get 


tr^Ct)  =  —  cr2(s)  + 


and  since  we  are  considering  the  time  when  s  =  o  this  gives 


a(t) 


s) 

V 


(20) 


Another  approach  to  this  problem  is  the  following.  At  each  station  i  we  may  deter¬ 
mine  the  apparent  time  of  fence  crossing  t  ^  which,  due  to  the  error  in  s . ,  will  have  an 
error  =  ct(s.)/v  as  above.  Then  the  best  estimate  t  can  be  found  by  our  standard 

method  and  will  have  an  error 


r  1 

1  1 

-  1 

r  1 

1  1 

1  ^ 

1 2) 

v2 

and  this  from  Eq.  (19)  gives  the  same  result,  a(t)  =  a(s)'v. 


Use  of  Phase  Rate 

A  proposed  development  in  the  Spasur  System  is  the  use  of  phase- rate  data  to  pre¬ 
dict  the  immediate  future  path  of  a  satellite  and  to  obtain  approximate  orbital  elements. 
We  begin  with  the  theoretical  procedure  for  such  calculations  and  then  discuss  practical 
simplifications  and  results.  Recall 


A.(p{rev)  =  rf  sin 
(rev/sec)  =  d  cos  fl  (rad/sec) 

from  which  we  obtain  the  error  formulas 

Velocity  Component  Determination  -  Velocity  components  from  angular  rate  are 
given  by 


Also 


The  most  accurate  values  of  <t>  and  <i>  are  taken  from  the  longest  baselines,  where 
Vd  is  small;  in  fact  V/d  -  lO'*.  This  means  that  except  for  very  low  or  very  fast  sat¬ 
ellites  we  may  ignore  the  second  term  in  the  expression  for  Similarly  we  may 


Fig.  6a  -  Level  contours  of  <r(a)/<T(T])  x  10'^  or  *  10'*:  stations  1  and  2.  The  curve  parameters  nuiy  be 

interpreted  as  error  in  statute  miles  or  statute  miles  per  second  assuming,  for  example,  data  having  18  degrees  (50 
coiints)  rms  phase  or  phase-rate  error  from  a  455-ft  (50  baseline. 


Fig.  6b  -  Level  contours  of  ofiVo’f’?)  »  10**  or  <r(vi)/‘^(v)  “  1®**:  stations  1  and  3,  The  curve  parameters  may  be 
interpreted  as  error  in  statute  miles  or  statute  miles  per  second  assuming,  for  example,  data  having  18  degrees  (50 
counts)  rms  phase  or  phase- rate  error  from  a  455- ft  (50  X)  baseline. 
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usually  ignore  the  second  term  in  the  expression  for  Thus  we  have  the  following 

basic  formulas; 


■  J. » 

■’<*>  =  c.b 

(21) 

=  r,  (r(0)  . 

North- South  Velocity  Component  -  As  was  true  in  the  analysis  using  phase  in  the 
subsection  "North-South  Errors,"  the  phase-rate  data  at  any  one  station  is  sufficient  to 
determine  velocity  once  the  position  "in  the  fence"  is  known.  Using  the  same  methods 
we  have  the  best  estimate 


V 


i 


-  1 


i 

z: 


i 


Using  )  =  'i  where  ^  is  the  north-south  space  angle  rate  we  have  for  two 
stations  ‘ 

(\  - 1  r  ^  r  ^ 

-7+-^)  (22) 

f,  T^j  r,  +  fj 

which,  as  stated,  is  the  same  form  as  the  north- south  position  error. 

In- Fence  Velocity  Component  -  With  the  data  from  two  stations  (Fig.  7a)  we  may 
determine  the  component  of  velocity  in  the  fence.  At  each  station  i  we  measure  the 
component  of  velocity  v.  perpendicular  to  r,,,  the  line  from  station  i  to  the  satellite 
(Fig.  7b).  The  problem  is  to  find  the  velocity  in  the  fence  Vf  which  has  these  compo¬ 
nents.  To  do  this  write  each  7.  as  a  magnitude  v. ,  which  we  measure,  times  a  unit 
direction  vector,  which  is  known  from  the  position  of  the  satellite. 

=  Vj  [- sin  (0j  +  a),  cos(0j  +  a)] 

V2  =  V2  [- sin  (^2  ”  cos(02”'^)J  ■ 


y 


(a)  Two -station  fix  (b)  Velocity  in  the  fence,  v, 

Fig.  7  -  Geometry  of  the  in- fence  velocity  determination 
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We  express  the  fact  that  vj  is  the  component  of  v,  in  the  direction  of  using  the  vec¬ 
tor  dot  product: 


Vf  •  [-sin(Hj  +  a),  cos(^'j  +  a)l  =  V| 


Vf  •  I-  sin  ff'j  -  a) ,  cosff'j-o  ]  =  V2  • 


Then  substituing  v,  =  [v^,  v^]  in  these  equations  we  get  two  simultaneous  linear  equa¬ 
tions  in  and  Vy  which  we  may  solve  to  find 

v_  =  — :  [v|  cos  (■(’,  -  a)  -  V,  cos  ((',  +  a)] 

*  sin  /5  *  ^  ^  * 


V 


y 


1 

sin  fi 


Iv, 


sin(('2  -  -  Vj  sin^Hj  +al] 


2 

Vf 


2  V j  V2  cos 


(23) 


where  -2a,  In  vector  form 


|vi  -  V2I 

I  sin  til 


From  Eq.  (23)  we  obtain  the  error  formula 


<a^(Vf)  =  -  y  - -  [(Vj  -  Vj  cos  C7^(vy)  +  (vj  -  Vj  cos  ti)^ 


Vj  sin^/^ 


+ 


I  V2 

5in2^\  ''f 


Vf  cos 


ct*(/5> 


where 


Total  Velocity  -  Finally  for  total  velocity  v  we  have 


(24) 


222 

V*  =  Vj  + 


V*  cr*(v)  =  V  *  a2(v^)  +  v/  (7*(v^). 


(25) 


These  components  must  be  corrected  for  the  effects  of  the  earth's  rotation,  but  for  the 
purpose  of  the  following  error  discussion  this  change  of  coordinates  may  be  ignored.  It 
is  now  possible  to  start  with  phase  and  phase- rate  information  from  two  stations  together 
with  the  standard  deviation  in  these  quantities  and  calculate  a  position  vector  and  a  veloc¬ 
ity  vector  and  also  the  standard  deviation  in  the  magnitude  of  these  vectors  or  their  com¬ 
ponents.  These  computations  are  somewhat  long;  we  will  consider  some  special  cases 
below. 


Orbital  Elements  -  In  our  discussion  of  orbital  elements  we  shall  use  the  following 
standard  units:  1  canonical  unit  of  distance  ^  1  earth  radius,  1  canonical  unit  of  time  s 
806.832  sec.  This  unit  of  time  is  chosen  so  that  if  g  is  the  gravitational  acceleration, 
then  gR^  -  1  (distance^  time  *2).  In  this  system,  formulas  for  conversion  from  position 
and  velocity  vectors  to  orbital  elements  do  not  involve  physical  constants.  If  we  let  a  be 
the  semimajor  axis  and  P  the  period,  we  have 
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r 


r  ^  2”  .t’  2 . 


(26) 


For  n  we  have  the  error  formula 


^  [  r)  +  (  ^  v)]  . 

^  4  ^  J 

(2-  rv2) 

We  will  assume  in  the  subsequent  examples  that  the  position  error  term  is  small  com¬ 
pared  with  the  velocity  error  term.  This  allows  us  to  write 


a)  =  2n^  V  ir(v)  . 


(27) 


Also  we  have  «(P')  =  3w  ^  „(n). 

In  theory  we  could  propagate  the  error  through  the  calculations  of  the  other  orbit 
parameters,  but  as  they  depend  on  the  direction  as  well  as  the  magnitude  of  the  velocity 
vector  they  are  correspondingly  more  difficult  to  compute.  A  straight  line  or  circular 
orbit  approximation  to  the  true  orbit  may  be  computed  with  only  these  elements  and  a 
third,  the  inclination  of  the  orbit  plane  to  the  fence  plane.  If  the  fence  plane  were  the 
equatorial  plane  this  would  be  the  parameter  usually  called  the  inclination  of  the  orbit 
plane.  The  error  introduced  by  this  approximation  in  the  near-future  predicted  position 
is  not  expected  to  be  larger  than  the  error  already  present.  Let  i  be  the  angle  between 
the  two  planes.  We  have  then  tan  i  =  from  which  we  derive 

<7(  i  )  =  (7^(Vf)  . 


Some  Test  Cases  -  We  will  consider  a  set  of  satellites  in  circular  orbits  whose  inter 
section  with  the  fence  occurs  midway  between  two  stations  10  degrees  apart.  Given  the 
orbital  radius  r  and  the  inclination  to  the  fence  i  a  circular  orbit  is  determined  and  we 
may  express  errors  in  the  orbital  elements  in  terms  of  phase  and  phase- rate  error.  See 
Table  1  and  Figs.  8  and  9.  Note  that  for  circular  orbits 


a) 

a 


cy(v) 

V 


and 


<7(P) 

P 


-  o^(v) 

•*  V 


The  last  term  in  Eq.  (24)  due  to  error  in  position  was  neglected  in  these  calculations. 


PRESENTATION  OF  RESULTS 

This  section  summarizes  the  results  of  our  calculations,  together  with  some  con¬ 
cluding  remarks.  The  error  functions  computed  are  defined  again  for  the  reader's  con¬ 
venience,  so  that  this  section  is  almost  self-contained.  Figures  3-5  consist  of  a  set  of 
equal  error  curves.  Each  curve  has  an  attached  number,  and  the  magnitude  of  the  error 
along  a  curve  is  equal  to  the  product  of  this  number  and  the  reciprocal  of  the  scaling 
factor. 


Explanation  of  Figs.  3,  4,  and  5. 

In  Figs.  3-5,  r  and  C  are  earth-centered,  earth-fixed  polar  coordinates  (see  Fig.  2). 
As  explained  in  the  preceding  section 


rms  error  in  r 


I 


28 


NAVAL  ftESEAftCH  LABORATORY 


Cfl 

u 


O 


u 

Si 


0 


u 


O 


«-*  c 

•04 


u 

o 


9 

01 


a: 


*■  ^ 
C  w> 

bS 

0 

o 

Oi 

eg^00^OCDe0O00C0^M^O9) 

^f-4^C^C0C0^lAiA(Dt«^00O>OO 

0 

o 

u 

••4 

^^coe^QomcMO)t*iAeA^09>oo 

o«H^e^esieo^^'iACDr^^oooiSo)0 

Sx' 

b 

o 

O 

O) 

f-4O>O)O>OiMco(OO>eor*CQr-esi0) 

^^MeAiAcof^ooo>*H<^^int*od 

^  ^  vH  ^  *H 

0 

o 

N 

•  •4 

cjtf4C0^i/><Dt*a0O<^4eg^'iAt*0) 
^  ^  ^  ^  ^ 

\  1 
^  tti 

.•4  0) 

'rs 

O 

S 

ooooooooooooooo 

ooeo>HO>(Dm^N(Oc<iO»Oe>30>^ 

mt~f-i04<^'Heor>«N4<coOr<- 

eoa>mr><-ito>-«(o>-ir~r^O>meMoo 

•-ii-iejegpsM^^irtirtcor-t- 

i 

ooooooooooooooo 

^eM(OomooNr-^Om>HiAt> 

<ooc-o^CMeM4<eOLn^^t-i-ic> 

^e->HiAomomO(Dc'iao^<-<e> 

0) 

(0 

.p 

'5' 

II  > 

0 

S 

^0>O0>t*^€D00lAe»Cb^C0^t» 

00CAOt*»tA^Ca^«H«H^CQC0^lA 

«-^c^)c^9eo^ir>(Dt*oO0^O^eMeo 

^  Vai4  ^ 

o 

O 

00 

^  (o  o>  r* 

00  CO  O  tn 

vH  eo  CM  CO 

o 

o 

CO  00  ^  A 

O)  ^  O  00  lA 
^  C^  CO  CO 

0 

o 

C4  t*  ^  lA  f-H 

O  ^  ^  00  CD 
^  ^  CM  CM  CO 

0 

o 

lA 

^  CM  lA  OO  ^ 

^  lA  ^  00  CD 
*H  OJ  CM  CO 

0 

O 

CD  00  D>  O  00 

CM  lA  ^  0>  CD 
^  CM  CM  CO 

0 

O 

to 

CO  •><  O  O  O 
m  (O  M  o  c> 

"H  »H  *<J  CCi  « 

o 

O 

M 

CO  oo  (O  CO  eg 
•4<coego>e- 
^  eH  eg  eg  eo 

p 

O 

00  O  00  t-  ^ 

^  t-  CM  0>  C» 

1-1  ^  CM  CM  eo 

i 

o<-io>ooir>e»eg4'coaoo>r>mr>co 

inr-ego>t-io^eocgegegco^ine- 

>H>-iegegcovmcDt>cooO'Hcgeo 

^  fM  ^  ^ 

H 

> 

b 

g 

b'5. 

OO^CD-^CplACOCM  -^oo>ooc-r- 
C-CMOOlACOOOOCD^CMOr-lACO-^ 
COCDaO-4^t*0>CMlAOO^eOCDO>eM 
^^^^C4CMeMCOCOCOCO^ 

«»■ 

> 

b 

hi 

OCDCDOOeMeMCD-^t-^eo^OO>0> 

OcDQO<OOC-^eMO>t-lACO-HQOCD 

t-r*o>eMiAr-ocoiAoo-H^t-o>eM 

-^-H^CMCMCMeMCOCOCOCO^ 

0^ 

’fm 

B 

eoaOT-im<Haot>cDCDe-o>i-<^rr-i-< 

co'^co^oeor-coin'ccocoeg.^r-i 

^^”4'^^eoeoe«seoeoeoeoeocoeo 

“•I 

egtootcoaoinegoooooaooi-Hcoco 

coooocoineo'M^iDOtcgmoegin 

•H^egcgeg'egc>»eocoeo«^'4<ir>io 

■“1 

r 

(earth 

radii) 

-^CMCO^O^lACDt-COOiOi-tCMeO^lA 

«H*H^4-4i-HiH^^^eMCMCMCMeMCM 

u 

0) 


fO 


♦Taking  a( i)  =  0.05  rev/sec,  \  =  278  cm,  and  d  =  1  km  we  have  cr(r))  =  (K/d)  a(ri>)  =  O.J39  *  10" 
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Fig.  8  -  Total  velocity  error  as  a  function  of  altitude  and 
inclination  to  the  fence 


Fig.  9  -  Error  in  inclination  to  the  fence  as  a  function 
of  altitude  and  inclination  to  the  fence 
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^  cT(i>)  y<ir*>  =  rtns  error  In  f,  assuming  no  error  in  position."" 

We  can  define  by  replacing  r  by  4  on  both  sides  of  these  equations.  Also, 

define  \/<Sw2>  by  the  equations 


we 


<7f(P) 


y<sw^>  =  y<s?  ■ 


^  cy( i)  y<Sw*>  =  \/<Sv  • 

Cl 


where  is  the  position  error  vector  (pointing  from  the  true  position  to  the  measured 
position),  and  is  the  velocity  error  vector  computed  after  assuming  no  error  in 
position. 

In  the  figures,  the  values  of  •J<hT‘^>  and  have  been  scaled  by  a  factor  of  10’^. 

For  a  value  of  k/a  of  2  «  10‘*  and  assuming  a(ip>  and  to  be  0.05,  corresponding  to 
5%  phase  noise,  we  have  values  of  (A./H)  a((p)  and  (Vd)  ad)  on  the  order  of  10'^.  Hence 
the  numbers  attached  to  these  curves  are  realistic  and  as  indicated  in  the  figure  titles, 
the  curve  parameters  may  be  interpreted  as  coming  from  phase  data  having  18  degrees 
(50  counts)  rms  phase  error  from  a  455-ft  (50  k)  baseline.  The  numbers  associated  with 
the  dots  on  the  vertical  axis  represent  the  minimum  (scaled)  values  of  the  errors. 

The  values  of  have  not  been  scaled,  so  that  the  actual  values  of  \/<S4^>  will 

never  be  more  than  a  few  hundredths  of  a  degree  in  the  region  covered  by  our  plots.  But 
although  these  angle  errors  may  be  small  numerically,  they  are  in  no  sense  "negligible," 
because  a  small  angle  error  in  an  earth- centered  polar  coordinate  system  may  be  blown 
up  many  times  when  the  earth- centered  system  is  transformed  to  a  coordinate  system 
centered  on  the  surface  of  the  earth.  In  other  words,  the  fact  that  is  small  does 

not  tell  us  anything  about  the  adequacy  of  Spasur  to  provide  search  coordinates  for  any  of 
our  tracking  systems,  since  none  of  our  tracking  apparatus  is  located  at  the  center  of  the 
earth. 


Explanation  of  Fig.  6 

The  same  set  of  curves  can  be  used  to  determine  both  north-south  position  and 
velocity  errors  and  are  found  in  Figs.  6a  and  6b.  Since  these  curves  are  also  scaled  by 
10'^  and  since  from  Eq.  (17),  <7(7?)  =  (X/d)  ct(ip)  and  ct(v)  =  (X/d)  cr((p),  the  values  on  these 
curves  are  also  realistic.  For  example,  a  point  on  the  level  curve  0.4  (statute  miles) 
represents: 

1.  An  error  in  north- south  position  of 

X  10-3  =  0.4  mile. 

a(r)) 


Therefore 


cr(s)  =  400  a(T)) 


or,  since  <7(7,)  »  10'^  radian. 


■i'We  remind  the  reader  that  in  the  calculations  we  set  (X,/d)  a(ifi)  -  (X/d)  (r(ii)  =  1.  The 
occurrence  of  y/<l>T^>  instead  of  y<Sf2>  in  the  second  equation  is  not  a  misprint.  See 
page  8  for  an  explanation. 
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c7(s)  =  400  *  10'®  -  0.4  mile. 
2.  An  error  in  north- south  velocity  of 


<y(v)  * 


10*^ 


0.4  mile. 


Therefore 


=  0.4  cr(T)) 

or,  since  cr(^)  *  10'^  radian/sec, 

=  400  X  10'^  =  0.4  mile/ sec. 


Explanation  of  Figs.  8  and  9  and  Table  1 

Properties  of  east-west  velocity  error  and  error  in  orbital  elements  have  been 
studied  for  the  special  case  of  satellites  in  circular  orbits  whose  intersection  with  the 
fence  plane  is  equidistant  from  the  two  stations.  We  have  taken  the  two  stations  to  be  10 
degrees  apart.  (Actual  earth- centered  angles  between  stations  are:  stations  1-2,  8.4 
degrees;  2-3,  13.4  degrees,  and  3-4,  8.1  degrees.)  In  this  situation  we  have  calculated 
the  percentage  error  in  total  velocity  for  satellites  at  various  altitudes  and  inclinations 
to  the  fence  as  well  as  the  error  in  the  measured  value  of  this  inclination. 

When  the  inclination  to  the  fence  i  =  0  degrees,  the  velocity  vector  lies  in  the  fence 
and  is  equal  to  the  "in  the  fence"  component  Vf .  When  i  »  90  degrees,  v  =  v^,  the  per¬ 
pendicular  component  discussed  above.  Note  in  Table  1  that  cr(v)/v  depends  very  little 
on  i .  This  is  probably  approximately  the  case  at  other  reasonable  positions  in  the  fence 
and  it  means  that  for  estimate  purposes  the  rather  messy  calculation  of  <7(v)  could  be 
replaced  by  ,  which  is  much  easier  to  compute.  Note  that  in  these  remarks  we  do 
not  consider  the  effect  of  length  of  time  in  the  fence;  cr(q>)  will  reflect  this  effect. 


Applications 

In  the  introduction  we  mentioned  the  use  of  error  as  a  weighting  factor  when  data  is 
received  from  more  than  two  stations.  The  phase- rate  information  on  velocity  may  be 
used  to  calculate  orbital  elements  after  a  single  pass.  Although  it  is  probably  unrealistic 
to  expect  these  values  to  be  very  good,  especially  in  the  case  of  the  more  involved  calcu¬ 
lations  such  as  eccentricity,  it  should  be  possible  to  obtain  an  estimate  of  period,  velocity, 
and  rather  good  estimates  of  inclination  and  right  ascension  which  may  be  used  to  deter¬ 
mine  a  rough  estimate  of  the  time  and  position  of  the  next  crossing.  Perhaps  a  more  use¬ 
ful  ai^lication  would  be  the  use  of  velocity,  inclination,  and  right  ascension  to  predict  the 
location  of  the  satellite  until  it  can  be  picked  up  by  a  tracking  radar.  For  this  purpose  we 
might  use  a  circular  orbit  approximation  (or  even  straight  line  path)  with  the  velocity  in 
this  assumed  orbit  taken  as  the  measured  v.  Thus  if  the  satellite  is  in  an  eccentric  orbit 
and  is  caught  at  perigee  the  velocity  will  be  greater  than  a  satellite  actually  in  the  circu¬ 
lar  orbit  assumed,  whereas  if  it  is  caught  at  apogee  the  velocity  would  be  less. 

Another  example  of  a  question  which  could  be  considered  in  the  light  of  the  error 
information  above  is  that  of  the  optimum  north-south  to  east-west  baseline  ratio.  Clearly 
if  cr(Vf)  is  much  greater  than  the  accuracy  of  the  latter  is  of  little  avail  in  com¬ 

puting  the  total  V.  We  should  like  these  two  quantities  to  be  roughly  the  same  magnitude. 
Specifically,  since 
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we  want 


_J _  ^  J 

V  (Vj^) 


In  the  special  case  of  satellites  intersecting  the  fence  midway  between  stations  10  degrees 
apart  this  gives 


^NS  ^ 

^EW  '^EN^ 


COS  cos 


2 

2 


where  is  the  velocity  component  measured  from  one  station,  and  where  vj^  v,  depends 
on  the  direction  of  passage  through  the  fence.  Since  the  satellite  is  in  the  fence  for  the 
same  length  of  time  viewed  by  north- south  network  as  by  east- west,  we  may  assume 

FroTn  these  considerations  we  would  conclude  that  Hf,s  and  should 
be  the  same.  If  the  same  argument  is  applied  to  position  error  it  turns  out  that 
should  be  longer.  For  instance  at  an  altitude  of  100  miles  midway  between  stations  1  and 

2  we  calculate  <T(r)  =  846fVd)<7^<p)  and  (t(s)  =  217(\'H)  ofep)  where  s  is  the  distance 
from  the  fence  plane  to  the  satellite.  At  200  miles  the  values  are  607  and  249,  and  at  250 
miles  and  100  miles  off  axis,  850  and  233.  Since  the  east-west  terms  are  on  the  order  of 

3  times  the  north-south  terms,  the  east-west  baseline  ought  to  be  about  3  times  the  north- 
south  to  get  compatible  values  for  these  particular  errors. 
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